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ABSTRACT 

We propose that the large, radius of ~ 10 3 AU, circumbinary rotating disks 
observed around some post-asymptotic giant branch (post-AGB) binary stars are 
formed from slow AGB wind material that is pushed back to the center of the 
nebula by wide jets. We perform 2D-axisymmetrical numerical simulations of fast 
and wide jets that interact with the previously ejected slow AGB wind. In each 
system there are two oppositely launched jets, but we use the symmetry of the 
problem and simulate only one jet. A large circularization-flow (vortex) is formed 
to the side of the jet which together with the thermal pressure of the shocked 
jet material accelerate cold slow-wind gas back to the center from distances of 
~ 10 3 — 10 4 AU. We find for the parameters we use that up to ~ 10~ 3 M Q is 
back-flowing to the center. We conjecture that the orbital angular momentum of 
the disk material results from the non-axisymmetric structure of jets launched by 
an orbiting companion. This conjecture will have to be tested with 3D numerical 
codes. 



1. INTRODUCTION 

Among the many riddles related to the transition from the asymptotic giant branch 
(AGB) to the planetary nebula (PN) phase, prominent is the presence of large circumbi- 
nary Keplerian disks around some post-AGB stars (Waters et al. 1993, 1998; Jura et al. 
1995, 1997, 2000; Van Winckel et al. 1995, 1998; Jura & Kahane 1999, Van Winckel 1999; 
Bujarrabal et al. 2003, 2005; Dominik et al. 2003; Hinkle et al. 2007). These post AGB 
stars (the primary) are found to have a close companion (the secondary) at a typical orbital 
separation of a ~ 1 AU, while the disk is much larger, up to ~ 2000 AU (de Ruyter et al. 
2006). In most cases the companion is thought to be a low mass main sequence star, but 
can also be a WD (e.g., Hinkle et al. 2007). In some systems the large-disk resides in the 
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center of a much larger bipolar nebula, e.g., the Red Rectangle (Cohen et al.2004) and 89 
Here (Bujarrabal et al. 2007). 

The total angular momentum in these disks is small compared to that of the binary 
system because the total mass in the disk is small < 1O~ 2 M (Bujarrabal et al. 2003, 2005). 
However, the specific angular momentum in the disk is about an order of magnitude larger 
than that in the binary system. This is not easy to explain (Soker 2000). One possible 
explanation is that the binary system exerts torque on the inner boundary of the disk and 
transfers angular momentum to the disk via viscosity (Jura et al. 2002; Frankowski 2007). 
It is not clear, however, if there is sufficient time for the viscosity to distribute angular 
momentum in such large disks (see section 3). This still leaves the question of how the 
material extends to distances of ~ 10 — 1000 times the orbital separation without dispersing 
first. After all, most of the gas is expelled from AGB and post-AGB stars. Motivated by 
the large specific angular momentum in these disks, Soker (2000) suggested that the dense 
matter in the equatorial plane is a dense slow equatorial flow rather than a Keplerian disk. 
However, the disk in the Red Rectangle was found to be Keplerian (e.g., Bujarrabal et al. 
2005), although it also has an outflowing component. 

In this paper we propose that the circumbinary disk is formed by the interaction of jets 
with the slow wind of the AGB or post-AGB primary star, that cause material to fall from 
large distance back toward the center. In this preliminary study we present the principles of 
the flow structure, and don't try to fit specific cases. 

2. NUMERICAL SIMULATIONS 

The simulations were performed using Virginia Hydro dynamics- 1 (VH-1), a high-resolution 
multidimensional astrophysical hydrodynamics code (Blondin et al. 1990; Stevens et al., 
1992; Blondin 1994). For gas temperatures above 10 4 K we use the radiative cooling func- 
tion for solar abundances from Sutherland & Dopita (1993; see Akashi et al. 2007), while 
for T < 10 4 K we use the radiative cooling time as given by Woitke et al. (1996; their fig. 
11). For numerical reasons we set a minimum temperature at 200 K. We use cylindrical 
(axisymmetrical) grid, namely, we simulate 3D flow with a 2D grid. Only one quarter of the 
meridional plane is simulated, as the other three quarters are symmetric to it. There are 
208 grids points in the 90° azimuthal direction and 208 grids points in the radial direction, 
with a cell size that increases with radius. 

We show here the results for two models. The first one has the following parameters. At 
the beginning of each simulation, t = 0, the grid is filled with the cold (1000 K) slow wind 
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material having an outflow velocity of v\ = 5 km s -1 . We assume that 950 years before the 
beginning of the jet-launching phase the slow wind mass loss rate increased by a factor of 
four. At t = we take the density for r < 3 x 10 16 cm to be that for a slow wind mass loss 
rate of M± = 2 x lO~ 5 Af yr _1 , and for r > 3 x 10 16 cm the density corresponds to a mass 
loss rate of Mi = 5 x 1O _6 M yr _1 . At t = we start to inject a conical jet within a half 
opening angle of a = 60° and a constant speed of Vj = 600 km s -1 . In recent years there 
are more indications for wide jets, in clusters of galaxies and in PNs (Sternberg et al. 2007; 
Sternberg & Soker 2008; Soker 2004), and therefore our use of a wide jet is justified. The 
inner boundary of the grid is at r = 10 15 cm. The jet is injected in a conical nozzle with its 
outlet at r = 7 x 10 15 cm. The mass injection rate of one jet is Mj = 10 _7 M Q yr _1 , injected 
uniformly within the angle a. In the second model we did not increase the mass loss rate of 
the slow wind prior to the jet launching. 

For numerical reasons we inject a very weak wind between the jet and the equatorial 
plane during the jet-launching phase. This outflow that must be incorporated to prevent 
numerical failure, prevents mass from being accreted directly to the center. For that reason, 
we examine the mass back-flow rate at larger distances. Therefore, the flow close to the 
center, where the equatorial slow outflow can be seen in the Figures, is not real. We rather 
expect an inflow there. 

In Figured] we show the entire nebula at t — 1900 years, where the colors represent the 
density and the arrows show the velocity. The corresponding gas temperature is shown in 
Figure [2J The structural features developed along the symmetry axis are not real. They may 
result form real physical instabilities, but their rate of growth and their exact structure are 
not real, but rather due to numerical limitations of the 2D numerical code. These numerical 
problems are limited to the vicinity of the symmetry axis. As we are interested in features 
near the equatorial plane, we will ignore the flow near the symmetry axis in this work, and 
postpone the discussion to the next paper. 

Most relevant to us is the flow of the 'cocoon' and the back-flow toward the center 
that is formed later close to the equatorial plane. The cocoon— a well known feature of 
expanding jets— is the slowly moving material around the expanding jet, which is formed 
from the post-shock jet material and some ambient matter (the slow-wind gas). In our 
simulations the cocoon forms a low density large circulating flow to the sides of the jets 
(in the axisymmetrical simulations the structure is a torus). The center of this vortex is 
seen at (x, z) = (1.6, 1) x 10 17 cm, where x is in the equatorial plane and z is along the 
symmetry axis. The jets and cocoon form two large low density bubbles, one on each side of 
the equatorial plane (in our simulation we only show one quarter of the meridional plane). 
As can be clearly seen there is a dense shell around each bubble. In a full 3D image we have 
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Density map and velocity field 




10 17 cm 

Fig. 1. — The density (scale on the right in log( g cm -3 )) and velocity plot in one quarter 
of the meridional plane at t — 1900 years after the onset of the jet, in the first model (with 
increased slow wind mass loss rate 950 years before jet launching). The arrows indicate 
the direction of the flow, and for three velocity ranges: v > 200 km s _1 (long arrow), 
20 < v < 200 km s _1 (medium arrows), and v < 20 km s -1 (short arrows). The vertical axis 
is the symmetry axis, and the horizonal axis is in the equatorial plane. 




Fig. 2. — The temperature (scale on the right in log( K)) map corresponding to the same 
run and time as in Figure [TJ 
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formed a bipolar nebula. 

As can be seen in Figures [T] and [2J the vortex touches a cool, T ~ 200 — 10 4 K, and 
dense material, seen as a stripe extending from (x,z) = (1.1,0) x 10 17 cm to (x,z) = 
(1.3,0.7) x 10 17 cm. Later this vortex (the circularization flow) and the thermal pressure 
of the hot bubble push cool material back to the center, as can be seen in Figure |3] that 
show the inner region of the grid at late times. This is the basic result of this and other 
simulations we performed with other parameters. Namely, that dense material can fall back 
toward the center from very large distances of fewxlO 17 cm. 

In the simulation presented here the mass of the back-flowing material increases with 
time. As stated, the total back-flowing mass is calculated by summing the inflowing gas 
mass at regions away from the slow equatorial outflow close to the center (this outflow was 
introduced for numerical reason). At t — 7000 years in the first model the total inflowing 
mass reaches a values of Mb ac k — 7 x 1O~ 4 M . For other parameters we get different masses, 
but of that order or smaller. 

The result of the second model (the one without increased slow wind mass loss rate) 
at t — 1300 year is shown in Figure HI The total (within the entire grid) back-flowing 
mass reaches values of Mback - 2-3 x 10 _4 M Q . The back-flowing mass within a radius of 
6 x 10 16 cm and 3 x 10 16 cm, is 2 x 1O~ 4 M and 2 x 1O~ 5 M , respectively. This is because 
at that time most of the back- flowing mass resides in the dense blobs seen in Figure H] in at 
distances of 4 x 10 16 cm <r<5x 10 16 cm. 

There are other features of the simulations that are less relevant to us. We mention 
two of these. As was proposed by Soker & Rappaport (2000), the pressure of the bubbles 
compresses material in the equatorial plane. The pressure is composed of both the thermal 
pressure of the relatively hot gas inside the bubbles and the ram pressure of the circularization 
flow (vortex). We emphasize this point: The jets, with the cocoon and bubbles, compress the 
dense equatorial gas. This is in contrast with many models which are based on preexisting 
equatorial dense gas that collimate the flow. The second interesting feature is the hot 
~ 1 — 3 x 10 6 K gas formed by the post-shock jet's material. This gas is a source of X-ray 
emission, which will be studied in a future paper. 

Precessing narrow jets will have the same effect as wide jets in inflating bubbles, if the 
precessing angle is large (~ 50°; Sternberg & Soker 2008). The same holds for narrow jets 
bent by the ram pressure of the slow wind, if they are bent by a large angle. 



- 7- 




Fig. 3. — The inner region of the simulated grid for two late times of the first model: 
4800 years (upper panel) and 5800 years (lower panel). Note the inflowing material. Shown 
are the density (scale on the top in g cm -3 ) and the velocity, represented by arrows for 
three velocity ranges: 10 < v < 100 km s _1 (long arrow), 5 < v < 10 km s -1 (medium 
arrows), and v < 5 km s _1 (short arrows). The vertical axis is the symmetry axis, and the 
horizonal axis is in the equatorial plane. Note that the dense equatorial outflow near the 
center, r < 2 x 10 16 cm, is introduced to prevent numerical failure. In reality we expect an 
inflow there, resulting from the back-flowing gas. 
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Density map and velocity field 




Fig. 4. — The density (scale on the right in g cm -3 ) and velocity plot in one quarter of the 
meridional plane at t — 1300 years after the onset of the jet in the second model (the model 
without increase in the slow wind mass loss rate). The arrows indicate the direction of the 
flow, and for three velocity ranges: v > 200 km s~ x (long arrow), 20 < v < 200 km s _1 
(medium arrows), and v < 20 km s~ x (short arrows). The vertical axis is the symmetry axis, 
and the horizonal axis is in the equatorial plane. Note that the dense equatorial outflow 
near the center, r < 2 x 10 16 cm, is introduced to prevent numerical failure. In reality we 
expect an inflow there, resulting from the back-flowing gas. 
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3. ANGULAR MOMENTUM 

We suggest that the back-flowing material discussed above forms the equatorial large 
disk found around close binary post-AGB stars. But where does the required angular mo- 
mentum come from? In the binary model for shaping planetary nebulae there are always 
departures from axisymmetry (Soker & Rappaport 2001). The jets are launched by the 
companion that is rotating around the center of mass. Two effects cause the jet's axis to be 
inclined to the z direction (symmetry axis). First is the orbital motion. Let us consider a 
jet launched by the secondary star at a speed of 400 — 800 km s _1 with its symmetry axis 
perpendicular to the equatorial plane in the secondary frame of reference. With a secondary 
orbital velocity of ~ 20 — 30 km s _1 the jet axis relative to the nebula will be deflected 
from the z axis by ~ 1 — 5°. In addition, the jet will be bent by the ram pressure of the 
primary stellar wind (Soker & Rappaport 2001). An eccentric orbit, as found in many of 
these systems, increases the effect of a temporal departure from axisymmetry. 

As a result of the bending, the inflated bubble will not be completely axisymmetric 
at any given time. In addition to the ram and thermal pressures toward the equatorial 
plane, there will be a small toroidal pressure gradient exerted by the bubble in the direction 
of orbital motion as the bent jets rotate around. This component will be very small, but 
nonetheless, sufficient to form a Keplerian disk because the material forming the disk in our 
model falls back from a large distance of thousands of AU to hundreds of AU, where the 
disk forms. The material escaping the system caries angular momentum with opposite sign, 
such that the total angular momentum in the system is conserved. 

Our conjecture that this departure from pure axisymmetry can lead to Keplerian disk 
formation should be examined with 3D numerical simulations, which is a topic of a future 
study. At this point we can only demonstrate the conjecture with the following estimate. 
Consider that the thermal pressure gradient and ram pressure accelerate matter to speed of 
v± < v a < Vj. We scale with v a = (viVj) 1 / 2 ~ 50 km s _1 , as we also find in the simualtions. 
The matter flowing toward the equatorial plane will be decelerated by the collision with 
matter coming from the other side of the equatorial plane, and its motion perpendicular 
to the orbital plane will vanish on average. On the other hand, the toroidal acceleration is 
in the same direction on both sides of the equatorial plane. Let the toroidal acceleration 
lead to a toroidal speed of v t = x v a at distance x t from the center in the equatorial plane. 
The bent jets sweep around at a 'phase-speed' of uox t) where u = 2n/P, and P ~ 1 year 
is the orbital period. We take x ~ Vj/ivxt, which for xt ~ 10 4 AU and Vj ~ 600 km s" 1 
gives x ~ 0.002, hence Vt ~ 0.002t> a ~ 0.1 km s _1 . The specific angular momentum of the 
compressed mass is then, for the above values, jt = VtXt ~ VjV a /ou ~ 10 3 km s" 1 AU. For 
a central object mass of 1M Q the corresponding Keplerian orbit is at x K ~ 10 3 AU. We 
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conclude that in principle a departure from temporarily axisymmetry, as expected in binary 
systems, can lead to non-zero angular momentum of the bound equatorial gas formed by a 
back-flowing material. 

If the bubbles are inflated by precessing jets, or jets bent by the ram pressure of the 
slow wind, then the specific angular momentum of the compressed torus will be even larger. 

The viscosity interaction time scale to spread material in geometrically thin disks is 
t d ~ 0.1R 2 /a d C s H (Frank et al. 1985), where R is the radius of the disk, H ~ (C s /v K )R 
is the vertical size of the disk, vk(R) is the Keplerian velocity at R, and a d ~ 0.01 — 1 is 
the viscosity parameter. For H/R < 0.3 and a central binary mass of ~ 1M , we find the 
viscosity time to be t d > 5000(i?/1000 AU) 3 / 2 /a year. This is not much shorter than the age 
of post-AGB stars. In our model, therefore, a disk is built, but it does not reach a complete 
equilibrium, and does not settle into a fully thin disk. We note that this time scale does 
not allow for angular momentum to be transport by viscosity from the binary system to the 
disk. 



4. SUMMARY 

We propose that the large circumbinary rotating disks found around post-AGB stars 
result from slow-wind material that is pushed back to the center by wide jets that interact 
with the slow AGB or post-AGB wind. A large circularization and back-flow, the 'cocoon', 
is formed, which together with the thermal pressure of the shocked gas compress equatorial 
gas. Performing 2D numerical simulations we found that the mass of the bound material 
can reach values of ~ 10~ 4 — 1O~ 3 M , and it flows-back from distances of ~ 10 3 — 10 4 AU. 

We emphasize that in this paper we don't try to explain all forms of bipolar PNs, and 
not all morphological features. In particular, we are not aiming at explaining nebulae with 
narrow lobes, which require narrow jets (Lee & Sahai 2003; Dennis et al. 2007) and for which 
the compression of the material in the equatorial plane might require a different approach 
(Sahai et al. 2005). We also did not try to explain expanding dense equatorial torii or small 
disks having sizes of < 100 AU for which the viscosity time scale is shorter than the age of 
the system. We are aiming only at explaining bound rotating material at ~ 10 3 AU from 
the central binary system. 

The main ingredients of the proposed model are as follows. 

1. A companion that accrete mass and blow wide jets. The jets cannot be too fast or too 
slow. Speeds of ~ 600 km s -1 , as appropriate for solar-like and somewhat less massive 
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main sequence stars, work the best. 

2. During the jet-launching phase, the companion must accrete a large fraction of the 
mass lost by the primary AGB or post-AGB star, in particular in the equatorial plane. 
Otherwise, the outward flowing wind will destroy the disk. This implies that the 
companion should be close to the envelope, but outside the envelope. Either the mass 
transfer occurs via a Roche lobe overflow, or the companion resides in the acceleration 
zone of the primary wind. 

3. For the jets to form a large 'cocoon' the jets cannot expand as narrow jets to large 
distances. This implies that the jets should have a large opening angle (wide jets), or 
be narrow and precess fast or be bent by the ram pressure of the slow wind very close 
to the center. 

4. Because of (i) the orbital motion, (ii) the primary stellar wind that hits the jets 
launched by the secondary star, and/or (Hi) eccentricity, the axes of the two jets 
are not exactly perpendicular to the orbital plane. The two jets are bent relative to 
the symmetry axis in the same direction. This direction changes periodically as the 
companion orbit the center of mass, and might lead to small component of toroidal 
acceleration of the mass flowing back to the center. With a crude analysis we found 
that it is enough that this toroidal acceleration component is ~ 0.001 — 0.01 times the 
magnitude of the acceleration toward the equatorial plane to form a large Keplerian 
disk. The specific angular momentum of the compressed equatorial matter will be large 
in cases where the bubbles are inflated by precessing or bent narrow jets. Future 3D 
numerical simulations will have to examine the conditions for the bound material to 
posses enough specific angular momentum to form a large Keplerian disk. 

5. The viscosity time scale of large disks (~ 1000 AU) is not much shorter than the age of 
these systems. We therefore don't expect large post-AGB disk to be in full equilibrium. 

There are several questions that our study could not answer, and must be addressed 
before our proposed flow structure can be accepted as a possible explanation for the large 
disks observed around most post-AGB binary stars with orbital period of ~ 1 year. Firstly, 
the conditions for the back-flow to occur must be quite common. Namely, these binary 
systems must interact in a way that will lead to wide jets with velocity of Vj ~ 400 — 
800 km s^ 1 . Secondly, the jet interaction with the slow wind forms a large nebula. We 
predict that a large nebula exists around each of these systems. However, in some cases the 
nebula will be very large, and therefore of low density and below detection limits. Thirdly, 
parameters appropriate for specific nebulae should be examined. For example, it should be 
demonstrated that bound mass of ~ 0.01M Q (Bujarrabal et al. 2003, 2005) can be obtained 
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to explain the circumbinary disk in the red rectangle; we obtained a mass of fewxlO~ 4 M Q 
for the parameters we used in several simulations. This requires an extensive study, as there 
are many unknown in the process of jet- launching. 

We thank John Blondin for his immense help with the numerical code and for his 
comments on the manuscript. This research was supported by the Asher Space Research 
Institute. 
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